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ABSTRACT 

The POINT- AGAPE survey is an optical search for gravitational microlensing events towards 
the Andromeda Galaxy (M31). As well as microlensing, the survey is sensitive to many dif- 
ferent classes of variable stars and transients. In our first paper of this series, we reported the 
detection of 20 Classical Novae (CNe) observed in Sloan r' and i' passbands. 

An analysis of the maximum magnitude versus rate of decline (MMRD) relationship in 
M31 is performed using the resulting POINT- AGAPE CN catalogue. Within the limits of the 
uncertainties of extinction internal to M3 1 , good fits are produced to the MMRD in two filters. 
The MMRD calibration is the first to be performed for Sloan r' and i' filters. However, we 
are unable to verify that novae have the same absolute magnitude 15 days after peak (the ti$ 
relationship), nor any similar relationship for either Sloan filter. 

The subsequent analysis of the automated pipeline has provided us with the most thor- 
ough knowledge of the completeness of a CN survey to-date. In addition, the large field of 
view of the survey has permitted us to probe the outburst rate well into the galactic disk, unlike 
previous CCD imaging surveys. Using this analysis we are able to probe the CN distribution 
of M31 and evaluate the global nova rate. Using models of the galactic surface brightness of 
M31, we show that the observed CN distribution consists of a separate bulge and disk popu- 
lation. We also show that the M3 1 bulge CN eruption rate per unit r' flux is more than five 
times greater than that of the disk. 

Through a combination of the completeness, M3 1 surface brightness model and our M3 1 
CN eruption model, we deduce a global M31 CN rate of 65±^ year 1 , a value much higher 
than found by previous surveys. Using the global rate, we derive a M31 bulge rate of 

year -1 and a disk rate of 27l]J year -1 . Given our understanding of the completeness and 
an analysis of other sources of error, we conclude that the true global nova rate of M3 1 is at 
least 50% higher than was previously thought and this has consequent implications for the 
presumed CN rate in the Milky Way. We deduce a Galactic bulge rate of 14_g year -1 , a disk 
rate of 20l^ year -1 and a global Galactic rate of 34^^ year -1 . 

Key words: novae, cataclysmic variables - galaxies: individual: M31 
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1 INTRODUCTION 

Classical novae (CNe) undergo unpredictable outbursts with a total 
energy that is surpassed only by gamma-ray bursts, supernovae and 
some luminous blue variables. However, CNe are far more com- 
monplace than these other phenomena (Warner 1989). 

Since their first recorded observations CNe have subsequently 
been identified as a sub-class of cataclysmic variables (CVs). The 
canonical model for CVs (Crawford & Kraft 1956) is a close binary 
system, containing a massive C-0 or O-Mg-Ne white dwarf (the 
primary) and a low-mass near-main-sequence late-type dwarf that 
fills its Roche lobe (the secondary). Any increase in size through 
evolutionary processes of the secondary results in a flow of ma- 
terial through the inner Lagrangian point into the primary's lobe. 
The high angular momentum of this transferred material causes it 
to form an accretion disc around the white dwarf, whilst viscous 
forces within the disc act to transfer the accreted material inwards, 
resulting in a small amount of the accreted hydrogen-rich material 
falling on to the white dwarf's surface (King 1989). In CN systems 
the mass accretion rate is generally lower than 1O _9 M0 year -1 
(Cassisi et al. 1998). As the accreted layer grows, the temperature 
at the base of the material increases. Hydrogen burning in the ac- 
creted envelope soon develops. Given the correct conditions, this 
can lead to a thermonuclear runaway (TNR) in which the accreted 
enveloped (and possibly some of the "dredged-up" white dwarf) is 
expelled from the system in a nova eruption (King 1989; Starrfield 
& Iliadis 2005). 

CNe typically exhibit outburst amplitudes of ~ 10 — 20 mag- 
nitudes and display an average absolute blue magnitude of Mb = 
—8 at maximum light, with a limit of around Mb = —9.5 for the 
very fastest (Shara 1981b; Warner 1989). The ability to accurately 
measure the distance to many Galactic novae (using expansion par- 
allax techniques) and a correlation between a nova's luminosity at 
maximum light and its rate of decline (Hubble 1929; McLaughlin 
1945) makes them potentially useful as primary distance indica- 
tors. However, until recently, generally poor light-curve coverage, 
small sample sizes and a current lack of understanding of how the 
properties of CNe vary between different stellar populations, have 
severely limited their usefulness as standard candles. Nevertheless, 
their relatively high frequency allows novae to be used as a tool for 
mapping the spatial distribution of the population of close binary 
systems in nearby galaxies. CNe may also be used to test nuclear 
reaction models and theories, whilst nucleosynthesis during a nova 
eruption is thought to make a substantial contribution to the abun- 
dances of a number of chemical species in the Galaxy such as 13 C, 
15 N and 17 (Jose 2002). 

The "speed class" of a CN is often used to describe the overall 
timescale of an eruption and to classify a nova (McLaughlin 1939; 
Bertaud 1948). The definition of the various classes depends on 
the time taken for a nova to diminish by two (or three) magnitudes 
below maximum light, ti (or i 3 ). Throughout this paper we will 
use the speed class definitions given in Warner (1989). 

1.1 Maximum Magnitude, Rate of Decline relationship 

From his years of observations of CNe in M31, Hubble (1929) 
noted that the brighter a nova appeared at maximum the more 
rapidly its visible light diminished. Given that all M3 1 novae can 
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be considered to lie at equal distance from the observer, Hub- 
ble's observation clearly implied a relationship between the nova 
speed class and its maximum magnitude. These observations for 
extragalactic novae were later confirmed for Galactic novae by 
McLaughlin (1945), who used a combination of expansion paral- 
lax, interstellar line strengths and Galactic rotation methods to mea- 
sure the distances of the nearby novae. Over time the empirically 
determined maximum magnitude versus rate of decline (MMRD) 
relationship for CNe has become accepted and refined (Pfau 1976; 
de Vaucouleurs 1978; Cohen 1985; Downes & Duerbeck 2000). 

A recent calibration of the MMRD relationship was made by 
Downes & Duerbeck (2000) using new distances, derived from ex- 
pansion parallaxes, for a sample of 28 Galactic novae, given by 

M v = (-11.32 ± 0.44) + (2.55 ± 0.32) log (t 2 /days) (1) 

Downes & Duerbeck (2000) concluded that a linear relationship 
is sufficient to model the Galactic MMRD. They also derived a 
typical scatter of ~ 0.6 magnitudes for CNe about their linear fits. 
Much of this scatter is thought to be due to difficulties in measuring 
accurate distances to the novae (Gill & O'Brien 2000; Warner 2005; 
Shafter 2005) and from intrinsic scatter in the optical decline due 
to variations in outburst parameters. 

However, it is also known that the linear MMRD relationship 
is not valid for the fastest and slowest novae (Arp 1956; Schmidt 
1957). Novae from M31 and the Large Magellanic Cloud are bet- 
ter described in terms of a "stretched" S-shaped curve. The form is 
somewhat supported by theoretical modelling of the nova eruption 
(Livio 1992). The "flattening" of the MMRD for brighter novae is 
thought to be caused as the mass of the white dwarf in the cen- 
tral system approaches the Chandrasekhar limit (Livio 1992). Con- 
versely, the flattening of the MMRD for the fainter novae is thought 
to be an observational selection effect (Warner 1995). 

Capaccioli et al. (1989) were drawn to the conclusion that 
the same MMRD relationship is valid in all galaxies of all Hubble 
types. This idea can be exploited to combine data from many differ- 
ent galaxies. As a result, the MMRD relationship can be used as a 
fundamental distance indicator (Shara 1981a). However the use of 
the MMRD relationship as a viable distance indicator is dependent 
upon being able to accurately measure the maximum brightness of 
a particular nova and its speed class, requiring good sampling of 
both the maximum light and the decline. 



1.2 Absolute Magnitude 15 days after peak 

Buscombe & de Vaucouleurs (1955) observed that all CNe ap- 
peared to reach approximately the same absolute magnitude 15 
days after their maximum light (Mis). The apparent constancy and 
value of Mis is yet to be fully explained, despite attempts to place 
it on a more physical footing (Shara 1981b). Sometimes referred to 
as the £15 relationship, a recent calibration was carried out by Fer- 
rarese et al. (2003) using nine newly discovered novae from Hubble 
Space Telescope (HST) observations of M49. Their calibration is 

Mis.v = -6.36 ± 0.19 (random) ± 0.10 (systematic) (2) 

However, there is great inconsistency in the calculated values of 
M15 (see Table 2.4 of Warner 2005). More recent results have 
called into question the reliability of using this so-called "iir," rela- 
tionship for distance derivations and the validity of the relationship 
itself (Jacoby et al. 1992; Ferrarese et al. 2003). 
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1.3 CN in M31 and the global nova rate 

A large number of CN surveys in M31 have been carried out, re- 
sulting in the discovery of around 500 novae. These have indicated 
the global nova rate in M31 to be ~ 30 - 40 yr _1 (Shafter & 
Irby 2001). Table 1 in Darnley et al. (2004, hereafter Paper 1) sum- 
marises the findings of many of these past surveys. The relatively 
high nova rate in M31 and its close proximity to our own Galaxy 
are major advantages of targeting M31 for nova surveys. However, 
since M3 1 is nearer edge-on than face-on, with an inclination angle 
of ~ 77° (de Vaucouleurs 1958), the task of unambiguously dis- 
tinguishing between novae erupting within the disk or within the 
bulge is rather difficult (Hatano et al. 1997). Consequently there re- 
mains debate surrounding the distribution and rate of novae within 
M31. The large inclination angle also introduces additional extinc- 
tion complications. 

The early M31 CNe surveys of Arp (1956) and Rosino (1964) 
found that the nova distribution decreased significantly towards the 
centre of the bulge, with Rosino (1973) reporting the centre of the 
bulge to be "devoid of novae"; all of this was despite their attempts 
to detect novae within the central bulge regions. However, the first 
M31 Ha survey (Ciardullo et al. 1987) found that the nova distri- 
bution follows the galactic light all the way into the centre of the 
bulge. A combination of the Arp (1956) novae with the Ciardullo 
et al. (1987) catalogue yielded the result that the bulge nova rate 
per unit B light was an order of magnitude greater than that of the 
disk, implying that the vast majority of the M31 novae arise from 
the bulge population. This result was later confirmed by Capacci- 
oli et al. (1989) after undertaking a comprehensive analysis of all 
M3 1 CN data. However, there is the potential for biases due to ex- 
tinction, especially within the disk, as the Ha surveys had focused 
primarily on the bulge, using much earlier B band surveys to "fill 
in" the disk data. In an attempt to tackle the lingering extinction 
issues Shafter & Irby (2001) extended the Ha observations into the 
M31 disk. Using M31's planetary nebula distribution for compari- 
son, they arrived at the conclusion that the M3 1 CN distribution is 
consistent with an association with the bulge. 



1.4 Nova populations 

The idea that CNe may arise from two distinct populations was first 
postulated by Duerbeck (1990). This was further explored by della 
Valle et al. (1992) who presented evidence that fast novae were con- 
centrated closer to the Galactic plane than slower novae. Additional 
spectroscopic data have revealed that there may exist two spectro- 
scopic classes of CNe, the Fe II and He/N novae (Williams 1992). 
It has been shown the the He/N novae tend to cluster close to the 
Galactic plane and that they tend to be brighter and faster than the 
Fe II type (della Valle & Livio 1998). 

Theoretical studies of CN outbursts (e.g. Shara et al. 1980; 
Shara 1981a; Prialnik et al. 1982; Livio 1992; Prialnik & Kovetz 
1995) have shown that the form of the outburst depends upon prop- 
erties such as the white dwarf's mass, accretion rate and luminosity. 
These white dwarf properties may vary with the underlying stel- 
lar population. These findings lend support to the idea that CNe in 
differing stellar populations may have distinctly different outburst 
properties. 

The surface gravity of a white dwarf increases with increas- 
ing white dwarf mass. This leads to a higher pressure at the base 
of the accreted envelope when the TNR begins, resulting in a more 
powerful outburst. It also follows that, as the pressure at the enve- 
lope base is greater for more massive white dwarfs, a lower mass 



of accreted material is required for the envelope to achieve the tem- 
perature and density required for a TNR to be initiated. Thus, the 
more massive white dwarfs are expected to have shorter recurrence 
times and to exhibit faster light-curve evolution. 

1.5 The POINT-AGAPE CNe catalogue 

In Paper 1 we presented an automated pipeline that used objective 
selection criteria to detect and classify CNe within a dataset with 
good temporal sampling. We reported 20 CNe erupting within M3 1 
over three seasons, detected using the pipeline. Nine of these CNe 
were caught during the final rise phase and all were well sampled 
in at least two colours. The excellent light-curve coverage allowed 
us to detect and classify CNe over almost the full range of speed 
classes. 

For the purposes of the POINT- AGAPE microlensing survey, 
the Wide Field Camera (WFC) on the Isaac Newton Telescope 
(INT), situated at La Palma, was used to regularly monitor two 
M31 fields between August 1999 and January 2002. The field cen- 
tres were located at a = h 44 m 00>!0, 8 = +41°34'00'.'0 and 
a = h 43 m 23!0, 8 = +40°58'15'.'0 (J2000) The WFC consists 
of a mosaic of four 2048 x 4100 pixel CCDs and the field locations 
are indicated in Figure 1 in Paper 1 . The field placements were pri- 
marily chosen to be sensitive to compact dark matter candidates, 
or Machos, which are predicted to be most evident towards the far 
side of the M31 disk (Kerins et al. 2001). The observations were 
conducted over three seasons in at least two broad-band Sloan-like 
pass-bands (usually r' and i', augmented with g' during the first 
season). The full distribution of observations can be seen in Table 3 
and Figure 2 in Paper 1. 

The outline of the current paper is as follows. In Section 2 we 
discuss our treatment of the internal extinction of M3 1 . Section 3 
presents our analysis of the MMRD relationship for M31. The tis 
analysis is presented in Section 4. In Section 5 we describe in de- 
tail our completeness analysis of the POINT-AGAPE dataset and 
the CN detection pipeline. In Section 6 we present our analysis of 
the CN population of M31. Section 7 details our evaluation of the 
global CN rate of M3 1 and finally we summarise and discuss our 
main findings in Section 8. 



2 EXTINCTION ESTIMATION 

In M31, as with most disk and spiral galaxies, the vast majority 
of the dust lies close to the disk plane and the extinction is patchy 
(Holwerda et al. 2005). As such, we expect that novae within M31 
should suffer a varying range of extinction, dependent upon their 
position in the plane of the galaxy and their line-of-sight displace- 
ment. These extinction uncertainties may be problematic when try- 
ing to analyse some of the global properties of the POINT- AGAPE 
CN catalogue, such as the MMRD and tu, relationships, as well as 
the completeness, CN population and nova rate. 

To compute the extinction across different parts of the M31 
disk we employ synthetic stellar models (Girardi & Salaris, 2001; 
Salaris, private communication) to estimate the extinction free in- 
tegrated colour (< r' — i >) of the disk. This is compared to the 
observed colour across the disk to allow us to compute an extinc- 
tion map. Assuming that the M31 reddening curve is similar to that 
in the Milky Way, the theoretical and observed colours yield an ex- 
tinction map for each band. The global average for the extinction 
maps are all re-scaled to give a global average extinction though 
the disk of A(i) — 0.8, corresponding to the typical value for Sb 
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galaxies found by (Holwerda et al. 2005, see their Figure 11) after 
transforming to Sloan magnitudes (Schlegel et al. 1998). Without 
this rescaling, our extinction values are systematically high, which 
we suspect may be because our synthetic models contain a larger 
fraction of bluer stars than is typically present in the M31 disk. 



3 THE MAXIMUM MAGNITUDE VERSUS RATE OF 
DECLINE RELATIONSHIP 

As noted earlier, the MMRD relationship is important as it can po- 
tentially be used as a tool to derive the distance to an extragalac- 
tic population of CNe by comparison to the MMRD in our own 
Galaxy. 

The CN decline rates originally estimated for the 20 detected 
POINT- AGAPE novae (see Table 6 of Paper 1) are re-computed for 
this analysis. To calculate the value of ti for each CN, we linearly 
interpolate between points on the decline of each lightcurve. This 
is carried out for both the r' and i' observations. The g' observa- 
tions are omitted as they are only available for seven of the detected 
novae. 

The measured maximum light and computed ti values for 
each CN candidate are shown in Table 1. The novae PACN-99- 
01, PACN-00-01 and PACN-01-01 are excluded because they are 
likely already to have been in decline at first observation; hence it 
is not possible to accurately determine the uncertainty in their max- 
imum light or decline rate. In addition, only a small portion of the 
light-curve of PACN-99-07 is sampled. As such, due to its erratic 
behaviour, we are doubtful whether the classification as a very fast 
nova is a true representation of this nova's speed class. Also, there 
are no i' data available around maximum light for PACN-99-05 
and PACN-99-06, so it is not possible to determine a ti value for 
these novae. Therefore, only 16 novae are used for the r' MMRD 
analysis and 14 novae for the i! analysis. 

An initial evaluation of the linear MMRD relationship for both 
the r' and i' data produces a poor fit, in the sense that the scatter 
of both distributions is much greater than that implied solely by the 
photometric uncertainties. 

3.1 Maximum light uncertainties 

In an attempt to further refine our MMRD relationship, and either 
reduce or help to explain the large scatter, we allow that the bright- 
est observation of each nova is only a lower limit to that nova's true 
maximum light. In the majority of cases each maximum observa- 
tion is straddled by observations from the following and preceding 
nights, leading to a small error in the assignment of the true max- 
imum light. We then estimate the maximum potential error on our 
measurement of the maximum light. Taking a good estimation of 
the general slope of each light-curve to be 2/t 2 magnitudes/day, we 
calculate the amount that each light-curve could have possibly in- 
creased in brightness between the two points straddling the bright- 
est observation. The maximum potential error induced by missing 
the maximum light for each CN is also shown in Table 1 . The large 
maximum light errors derived for PACN-99-07 and PACN-01-05 
are due to a combination of a fast decline rate and poor sampling 
around maximum light. In particular, the speed class assignment of 
PACN-99-07 is known to be suspect; a much slower decline time 
seems more fitting. As previously mentioned, this CN is excluded 
from this analysis due to the uncertainty of its decline rate. 

The photometric and maximum light uncertainties are com- 
bined by simple addition as the maximum light error is a sys- 



tematic rather than a random error. The actual maximum light is 
equally likely to lie at any point between the observed value and 
the computed potential maximum value. For the purpose of fitting 
the MMRD relationship, we "re-sample" the "observation" of the 
maximum light to be the mid-point of the computed range for each 
CN. 

We use a minimum absolute deviation method to fit the data, 
as our errors are no longer Gaussian and are dominated by uni- 
formly distributed systematic uncertainties. The scatter on these 
MMRD fits remains at ~ 0.6 and ~ 0.7 magnitudes for the r' and 
i' data respectively. This implies that the scatter in the MMRD rela- 
tionships is not significantly affected by missing the maximum by a 
day or two, as is the typical interval expected in the POINT- AGAPE 
data. The similarity between these fits and those performed with- 
out considering the maximum light uncertainties is also indicative 
of the apparent minimal effect induced by considering the possibil- 
ity of missing the maximum light of a CN. 

3.2 Extinction corrections 

Each CN's light-curve may still be affected by extinction within 
M3 1 and our own Galaxy. The Galactic extinction in the direction 
of M3 1 is well defined and relatively small compared to the maxi- 
mum potential extinction experienced within M3 1 . The extinction 
experienced by each nova's light depends upon the column density 
between the nova and the observer, with the maximum potential 
extinction dependent upon the CN's position. 

We use the extinction maps (see Section 2) to provide an esti- 
mate of the maximum potential extinction experienced by each of 
the 16 CNe for which valid r' decline rates could be computed and 
by each of the 14 novae with valid i' decline rates. The computed 
average extinctions are shown in Table 1. Figure 1 shows the distri- 
bution of extinction-corrected and maximum-magnitude corrected 
r' and i' data respectively. The base of the error bars represents 
the true maximum observed light, with the length of the bars rep- 
resenting the absolute range of the actual maximum light. As with 
the maximum magnitude errors, we assume that the extinction is 
equally likely to lie at any value between zero and the maximum 
estimate. Again, for this analysis, the three error sources are com- 
bined by simple addition. As the maximum light error and the ex- 
tinction error are both independent absolute maximum errors, the 
maximum error that can be experienced due to both of these sources 
is simply the sum of the two - in the direction of increasing lumi- 
nosity. We again assume that the best-guess maximum-light flux 
occurs equally distant between the observation and the extreme 
maximum error for the purpose of the MMRD fitting. The fits to 
the maximum-light estimates are 

m r , = (14.5 ± 1.3) + (1.5 ± 0.8) log (t 2 /days) (3) 

rrv = (14.5 ± 1.0) + (1.5 ± 0.6) log (i 2 /days) (4) 

The scatter in the final MMRD fits is ~ 0.7 and ~ 0.8 magnitudes 
for the r' and i' datasets respectively, comparable to the mean error 
size. The MMRD data are also analysed in the linear region of the 
"S-shaped" curve (0.5 < log[200d/t 2 ] «S 1.5), yielding 

rrv = (13.0 ± 2.2) + (2.5 ± 1.4) log (t 2 /days) (5) 

771;/ = (11.0 ± 2.4) + (3.9 ± 1.6) log (i 2 /days) (6) 

with a scatter of ~ 0.7 magnitudes about the fits for both the r' and 
i' data, which is again comparable to the mean error. 

All four of the MMRD relationships calculated (Equations 3- 
6) are shown in Figure 1, along with a recent Galactic calibration of 
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Table 1. r' and i' maximum observed magnitudes and corresponding times, maximum magnitude uncertainties and average extinction correction 
(in magnitudes) for each CN detected in the POINT-AGAPE data. 



Nova r '(*o ) *2(r') Estimated max Estimated *'(io) *2(*') Estimated max Estimated 

error on r' average r' error on i' average i' 

maximum light extinction maximum light extinction 



PACN-99-01 


16. 


53 


± 0.03 a 


30.50 






-0, 


67 


16. 


.39 


± 





.03" 


37.53 






-0 


.51 


PACN-99-02 


18. 


,91 


± 0.03 


99.48 


-0. 


03 


-0, 


.57 


19 


19 


± 





04 


58.09 d 


-0 


.04 


-0 


.43 


PACN-99-03 


17. 


79 


± 0.02 


59.62 


-0. 


02 


-0. 


60 


17, 


60 


± 


0. 


,04 


34.16 


-0 


06 


-0 


.46 


PACN-99-04 


18. 


41 


± 0.04 


164.39 d 


-0. 


02 


-0, 


52 


18 


.34 


± 





07 


87.25 d 


-0 


.02 


-0 


.39 


PACN-99-05 


17. 


70 


±0.04 


25.82 


-0. 


04 


-0. 


66 


















-0 


.50 


PACN-99-06 


16. 


,17 


±0.01 


20.30 


-0 


13 


-0. 


67 


















-0 


.51 


PACN-99-07 


18, 


1 ±0.1 


9.80 


-4. 


.2 


-0 


65 


18, 


02 


± 


0. 


04 


2.28 


-0 


88 


-0 


.49 


PACN-00-01 


17. 


73 


± 0.04 6 


38.65 






-0 


58 


17 


.58 


± 





,08 b 


11.88 d 






-0 


.44 


PACN-00-02 


18. 


15 


± 0.03 


198.55 


-0. 


01 


-0 


65 


18, 


.85 


± 





,05 


817.19 d 


-0 


()0 a 


-0 


,49 


PACN-00-03 


18. 


54 


± 0.03 


33.02 


-0. 


.06 


-0. 


67 


18, 


19 


± 


0. 


,04 


22.44 d 


-0 


.09 


-0 


.51 


PACN-00-04 


17. 


.61 


± 0.03 


30.65 


-0 


07 


-0 


66 


17, 


33 


± 


0. 


04 


36.44 d 


-0. 


06 


-0 


.50 


PACN-00-05 


17. 


30 


± 0.01 


59.21 


-0 


18 


-0 


65 


17 


.11 


± 


0. 


,01 


198.32 


-0 


06 


-0 


,49 


PACN-00-06 


17. 


09 


± 0.01 


13.85 


-0 


09 


-0. 


65 


16 


.64 


± 


0. 


.01 


13.44 


-0. 


14 


-0 


,49 


PACN-00-07 


19. 


53 


± 0.04 


55.21 


-0 


03 


-0 


48 


19. 


.48 


± 


0. 


05 


100.27 d 


-0. 


02 


-0 


.37 


PACN-01-01 


18, 


45 


± 0.02 c 


213. 12 d 






-0. 


61 


18, 


16 


± 


0. 


04 c 


330.86 d 






-0 


.47 


PACN-01-02 


17. 


14 


± 0.03 


22.06 


-0 


.10 


-0 


64 


16. 


,71 


± 





.04 


17.08 


-0 


12 


-0 


.48 


PACN-01-03 


17, 


30 


± 0.04 


143.71 


-0 


02 


-0 


.62 


16. 


,88 


± 





06 


66.21 d 


-0 


03 


-0 


.47 


PACN-01-04 


17, 


90 


± 0.03 


47.29 


-0. 


05 


-0. 


65 


17, 


38 


± 


0. 


,04 


37.24 


-0. 


06 


-0 


.49 


PACN-01-05 


15 


90 


± 0.01 


28.15 


-0. 


93 


-0. 


57 


15, 


61 


± 


0. 


01 


15.83 


-1 


68 


-0 


.43 


PACN-01-06 


17, 


38 


± 0.01 




-0 


.12 


-0. 


63 


16 


.88 


± 





.03 


38.24 d 


-0. 


.10 


-0 


.48 



a the light-curve of PACN-99-01 was visible at, or shortly after, maximum light in the first observational epoch of the first season, so it was only 
possible to put a lower limit on its maximum light. b as PACN-99-01 for the second season. c as PACN-99-01 for the third season. d it was not 
possible to follow these light-curves through two magnitudes below their observed maxima, so the value of <2 has been estimated from the general 
trend of these light-curves. 



the MMRD relationship (Equation 1 in Downes & Duerbeck 2000). 
The Downes & Duerbeck (2000) calibration has been translated to 
the distance of M31 using a distance modulus of 24.3 magnitudes 
(Welch et al. 1986), yielding 

m v = (12.98 ± 0.44) + (2.55 ± 0.32) log (t 2 /days). (7) 

It is clear that neither the r' or i' MMRD relationships provide 
enough data points to investigate in greater detail the true form (lin- 
ear of S-shaped) of the M31 MMRD relationship. As previous work 
(see Section 1.1) has specifically identified a linear region within 
M31's MMRD relationship, we will use the relationships defined 
within the linear region (shown in Equations 5 and 6) as the r' and 
i' MMRDs for M3 1 . 

We note that extinction is a rather weak factor in determin- 
ing the slope of the MMRD in the linear regime, where all the no- 
vae are intrinsically very bright [M(r')ltal7]. Indeed our MMRD 
slope determination with or without the extinction correction is es- 
sentially the same. 



4 THE tis RELATIONSHIP 

As discussed earlier, the ti$ relationship may also be useful in cal- 
culating the distance to a CN population. The majority of previous 
tis calibrations have been carried out using V-band data. How- 
ever, due to the restrictions of the POINT-AGAPE catalogue, we 
can only attempt calibration using r' and i' data. 

Using these data, an initial calibration of the t\;, relationship 
for the POINT- AGAPE catalogue indicates that the photometric er- 
rors alone do not account for the extent of the scatter of luminosities 
at 15 days following maximum light. 



As was attempted for the MMRD data, we can try to decrease 
or at least explain the scatter in these data by taking into account 
the line-of-sight extinction. Using the data in Table 1 to recalibrate 
each light-curve, we therefore reassess both the r' and i' tis rela- 
tionships for the POINT-AGAPE novae catalogue: 

m 15 y = 18.0 ±0.9 (8) 

m 15 y = 18.0 ± 1.0. (9) 

Figure 2 shows superpositions of the 16 r' and 14 i' re-calibrated 
light-curves that have well defined maximum light magnitudes. The 
light-curves are all plotted in units of time since maximum light. 
Each plot shows the light-curve behaviour for the first 50 days fol- 
lowing each eruption. There is clearly little or no convergence of 
these light-curves at time around 15 days. 

The inclusion of extinction and maximum light uncertainties 
has actually slightly increased the scatter of the ii 5 values. This 
implies that the scatter in luminosities 15 days after peak is not 
solely due to uncertainties in the luminosity of the nova, indicating 
that any tis relationship may not be valid for these pass bands. 



4.1 Comparison with previous results 

By assuming a distance modulus for M31 of 24.3 magnitudes 
(Welch et al. 1986), our computed tis values (see Equations 8 and 
9) might be compared with those found previously (see Table 2.4 
of Warner 2005): 

M 15 y = -6.3 ± 0.9 (10) 
M 15 v = -6.3 ± 1.0. (11) 
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log(200d/t 2 ) log(200d/t 2 ) 

Figure 1. The relationship between the r' brightness at maximum light and the decay rate (V2 (r') = log[200(i/<2 (»*')]) °f me 16 POINT- AGAPE CNe with 
well defined maximum lights and decay rates (left pane). Likewise for i' maximum light (right pane). The range of the error bars represents the maximum 
estimated range of the combination of the extinction, maximum light and observational uncertainties. The dashed line indicates an unweighted fit performed 
on all the data (Equation 3), while the solid line shows an unweighted fit performed on the data in the range 0.5 ^ log[200c!/t2 (?*')] ^ 1.5 (see Equation 5). 
The vertical grey line represents the "slow" boundary of the linear region of the MMRD (log[200d/t2 (»*')] = 0.5). The grey shaded region represents the 
best fit Galactic "S-shaped" MMRD, and the black dotted line shows the best fit Galactic linear MMRD - both these Galactic MMRD are derived for V data 
and have been transformed to the M31 distance. 




t (Days past maximum) t (Days past maximum) 

Figure 2. A superposition of the recalibrated r' light-curves of 16 of the POINT- AGAPE CNe (left pane) and the superposition of the recalibrated i' light- 
curves of 14 of the POINT-AGAPE CNe (right pane). The light-curves have been time-shifted so that the times of their observed maximum light are coincident. 
Each line represents the linear interpolation of the light-curves between observations. 



However, direct comparison between our results and previous re- 
sults can not be easily made. All calibrations of the £i 5 relationship 
to-date have been in "blue" bands, whereas our calibration is per- 
formed in "red" bands. Nevertheless, given that CNe are expected 
to become bluer with time, the result that our £15 luminosities are 
fainter than all but the most recent of the previous "blue" calibra- 



tions (see Table 2.4 of Warner 2005) is not surprising. The scatter in 
our results is very large, and larger than those found in previous sur- 
veys. For instance, in a recent HST study of M49 (Ferrarese et al. 
2003) the V-band £15 relationship was found to have a — 0.43 
magnitudes. It should also be noted that, as an elliptical galaxy, 
M49 does not suffer from large extinction problems. However, the 
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Figure 3. Plot of the distribution of r' magnitude scatter (black line) and 
the i' magnitude scatter (grey line) between observed nova magnitudes for 
a range of times following maximum light. 

large degree of the scatter in the POINT-AGAPE data can not be 
solely explained by the combination of maximum light and extinc- 
tion uncertainties, whose mean error is ~ 0.7 magnitudes for both 
the r' and i' data. 

As a final test, if the tw relation is valid, then we would expect 
a minimum in the scatter of the light-curves at or around 15 days 
after maximum light. Figure 3 shows a plot of the scatter between 
the POINT- AGAPE light-curves over a large range of time follow- 
ing maximum light for the r' and i' data. From inspection of this 
plot, it is quite clear that the light-curves of the POINT-AGAPE 
sample do not exhibit behaviour consistent with the existence of 
a ti5 relationship. However, the r' scatter does seem to exhibit a 
minimum at ~ 30 days after maximum light and the i' scatter is 
minimised ~ 35 days following maximum. However, these min- 
ima are coincident with the end of data sampling for a number of 
the light-curves, so may just be indicative of the temporal coverage 
of the POINT-AGAPE light-curves themselves. 



5 CN DETECTION PIPELINE COMPLETENESS 

Because we have have selected our nova candidates using objec- 
tive selection criteria, we can assess the efficiency of our selec- 
tion pipeline. This allows us to compute the completeness of the 
POINT- AGAPE CN catalogue and aids us in probing the underly- 
ing CN distribution and to compute a robust estimate of the global 
nova rate. To measure the completeness of the catalogue, we seed 
the raw POINT-AGAPE data with re-sampled light-curves of our 
20 detected CNe. We then re-run the entire CN detection pipeline 
on these seeded data to allow us to compute the proportion of re- 
covered light-curves. 

5.1 Creating test light-curves 

The seeded light-curves are positioned on a grid within the aligned 
image data stack (see Section 3 in Paper 1), with a grid spacing of 



15 pixels (5 arcmin). This grid spacing is chosen to allow the closest 
possible spacing of seeded objects, whilst minimising overlap of 
each star's point spread function (PSF). Each light-curve is seeded 
at a random eruption epoch, such that at least one point of the light- 
curve occurred between the first observational epoch and the final 
epoch. 

In order to seed the detected novae at any random time we 
linearly interpolate their light-curves between successive observa- 
tions. Whilst this works well when the timescale between observa- 
tions is small, it becomes less reliable when the gaps are larger. The 
largest gaps in the observations are usually of order 2 weeks, but in 
a few cases, when light-curves are followed across two seasons, 
these are up to six months. The two light-curves affected by this 
are PACN-00-02 and PACN-00-05 (see Section 5 of Paper 1), these 
light-curves being linearly interpolated across the seasonal gaps. 
Given the form of the light-curve of PACN-00-02, we expect this 
method to be relatively reliable as an estimate of the flux. However, 
given the predicted transition phase minimum for PACN-00-05, this 
estimate is much less reliable. 

5.2 Seeding the raw POINT-AGAPE data 

The generated light-curves are added to both the raw data and the 
PSF-matched data (see Section 3.2 of Paper 1) using the NOAO 
IRAF package environment 1 mkob jects, which scales the rele- 
vant image's PSF profile to the correct luminosity, then adds the 
scaled PSF profile to the data, recalculates the Poisson errors and 
combines these with the data. 

5.3 Re-running the nova detection pipeline 

The CN detection pipeline is re-run on the seeded data; however a 
number of stages of the pipeline are not used. As both the raw and 
the PSF-matched data are seeded independently 2 , the image align- 
ment, trimming, PSF-matching and background estimate stages are 
not required. 

The seeded PSF-matched data are run though the aperture 
photometry pipeline (see Sections 4.2 and 4.3 of Paper 1) to pro- 
duce the preliminary list of recovered light-curves. PSF-fitting pho- 
tometry is then performed at the position of each of the seeded 
novae recovered from the aperture photometry stage. These nova 
light-curves are then passed back through the "peak detection" 
stage of the pipeline. However, all of the pipeline stages that are 
related to the colour light-curves are ignored. We are able to ignore 
the colour criteria as these are solely introduced to distinguish CN 
light-curves from the light-curves of other objects that may have 
passed through the previous stages of the pipeline 3 . However, as 
we know that all 20 of the seeded light-curves are those of CN dis- 
covered in the POINT- AGAPE data, they have already passed the 
colour criteria. 

1 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc. under cooperative agreement with the National Science Foun- 
dation. 

2 In order to maintain consistency between the seedings in the raw and 
PSF-matched data, the same random seed is used to regenerate the Poisson 
noise for both CN seedings. 

3 The make-up of the POINT-AGAPE observation strategy makes it im- 
possible to seed i' and g' band data across observing seasons as there is 
minimal i' band data available for the 1st season and no g' band data avail- 
able for the 2nd or 3rd. 
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A seeded CN light-curve may "fail" the pipeline for any of the 
following reasons: 

(i) The object has been seeded at a location in the M31 field 
where, due to the brightness of the background and/or surrounding 
objects, it is impossible to make a 10a detection of the object at 
any epoch. 

(ii) There are not five consecutive detections, either because the 
object is too faint to detect or because it has been seeded such that 
there are not five observations in which the nova was visible. 

(iii) The observed "peak" of the seeded nova is not significant 
enough to pass the primary peak test. This is either due to the galac- 
tic background or because the nova has been seeded "low down" in 
its light-curve, i.e. the actual peak has not been seeded in any of the 
observations. 

(iv) A seeded CN light-curve can fail the periodicity test, the 
secondary peak height test or the "< 90% of points in peaks" test, 
if the nova has been seeded close to a region of the image that also 
varies significantly with time. This may be due to a nearby variable 
star, a region of bad pixels or a saturated object. Again, it is possible 
for a light-curve to fail this test if data around the actual peak of the 
nova has not been seeded. 

The numerical results of the completeness run of the CN 
pipeline are shown in Table 2. 

5.4 Completeness distribution 

In order to compute the completeness, we subdivide each CCD into 
1 arcmin grid squares, with each grid square containing 144 novae 
seed points. Depending upon the size of the trimmed CCDs (see 
Section 3.1 of Paper 1), the CCDs contain between 200 and 242 
grid squares. 

We first compute the completeness distribution of the CN 
pipeline; this distribution tells us the probability of the pipeline 
detecting a CN - of a type originally detected by the pipeline - 
at any position in the POINT-AGAPE fields, given that the CN in 
question is "visible" at least once between (and including) the first 
and last observations. This completeness calculation takes account 
of all of the major factors that affect the completeness of the de- 
tection pipeline, including the temporal distribution of the POINT- 
AGAPE observations, the galactic surface brightness, the variety of 
CN light-curve forms and any interference by foreground objects. 

The generated completeness map is shown in Figure 4. This 
illustrates that the completeness is relatively flat across both fields, 
within the noise, at a value between about 30% — 40%. However 
the completeness does decline towards the centre of the galaxy, as 
the galactic background begins to increase significantly. 



6 M31 CN POPULATION 

If we adopt the simplest assumption, that the nova distribution in 
M31 follows the light distribution (Ciardullo et al. 1987) then the 
probability of a CN erupting at a particular point within M3 1 is pro- 
portional to the flux at that position. By requiring that a given CN 
erupts within one of the two POINT-AGAPE fields, we can com- 
pute the probability of a CN erupting at a particular point within 
M31. However, from our completeness calculations, we also know 
the probability of detecting a CN in each grid square, given that a 
nova erupts within that square. Hence we can use this to calculate 
the probability of detecting a CN within the POINT- AGAPE fields, 




-20 20 

x (arcmins) 

i ^ : : : : : : : : ~h 

0.05 0.1 0.15 0.3 0.35 0.3 0.35 0.4 0.45 

Completeness 



Figure 4. The POINT-AGAPE CNe detection pipeline completeness distri- 
bution. Each rectangle represents one of the four INT WFC CCDs and the 
white circles indicate the positions of the 20 detected novae. The origin is 
the centre of M31 at (J2000) a = h 42 m 44?324, <5 = +41° 16' 08'.' 53 
(Crane et al. 1992). Also indicated are ten representative M31 "isophotes" 
from the surface photometry of de Vaucouleurs (1958), along with repre- 
sentations of the positions and sizes of M32 (within the southern field) and 
NGC205. 

given that a nova erupts within them. This detection probability is 
given by: 



where Pi is the probability of a CN erupting in a particular grid 
square containing a flux fi and e, is the computed pipeline effi- 
ciency in the grid square in question. Because of the uniformity of 
e, the distribution of detection probabilities closely resembles the 
galactic light. However, as the completeness drops slightly towards 
the galactic centre, this distribution has a slightly weaker central 
dependence than the flux. 

If M31 consisted solely of a single population of stars, then the 
"nova follows the light" distribution would be a good model of the 
CN distribution. However, with recent evidence pointing towards 
separate bulge and disk populations of novae (see Section 1.4), it is 
likely that the detection probability model (Equation 12) needs to 
be modified. 

6.1 Modelling M31's galactic light 

In order to calculate the overall CN rate of M3 1 or to investigate 
the possibility of separate bulge and disk CN populations, we need 
to be able to compute both the bulge and disk component of the 
light at any given point within the galaxy. In order to do this, we 
need to create a model of the flux distribution of M3 1 . To perform 
this modelling we subdivide each CCD using the one arcmin square 
grid system that we employed for the completeness calculations. 
In order to try to fit the disk or bulge components of the M3 1 
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Table 2. The effect of each stage of our selection pipeline upon the synthetic classical nova catalogue. These steps are described 
in Section 4 of Paper 1. 



Pipeline 
stage 


CCD1 


North-field 
CCD2 CCD3 


CCD4 


CCD1 


South-field 
CCD2 CCD3 


CCD4 


All 

CCDs 


>.ppHpH f^hipctc 
JCCliCli OUJCCla 


35,239 


35,376 


35,244 


35,376 


35,910 


35,239 


35,378 


34 864 


78? 626 


Ohippft: cppHpH within Hntn 

V_/LHGk,lL> W1L11111 Llillil 


18,291 


18,886 


18,314 


18,401 


14,185 


18,082 


18,861 


18 527 


143 547 


1 f\rr nhipntc 
1UO UUJCClS 


16,519 


17,569 


17,339 


17,107 


13,352 


16,216 


17,667 


17 049 


132 818 




Pipel 


ine 1st 


pass — aperture photometry 










5 consecutive detections 


13,457 


14,982 


14,945 


14,428 


12,451 


14,573 


15,849 


15,303 


115,988 


^> l piniidiy pedK 


13,078 


14,768 


14,660 


14,106 


12,221 


14,196 


15,409 


14 7SQ 


113 227 


PprinHipit^ tpct 

rciiumciiy icm 


12,979 


14,676 


14,550 


13,965 


12,156 


14,118 


15,191 


14 583 


112 218 


Primary peak height 


12,859 


14,564 


14,477 


13,784 


12,113 


13,945 


14,938 


14,391 


111,071 


Secondary peak height 


11,011 


11,642 


12,460 


11,253 


10,373 


11,926 


11,062 


11,179 


90,906 




Pipeline 2nd pass - PSF-fitting photometry 










5 consecutive detections 


9,981 


11,549 


12,330 


11,191 


10,157 


11,502 


10,628 


10,787 


88,125 


1 primary peak 


9,083 


11,056 


11,907 


10,696 


9,728 


10,837 


9,978 


10,157 


83,442 


Periodicity test 


9,083 


11,056 


11,907 


10,696 


9,728 


10,837 


9,978 


10,157 


83,442 


Primary peak height 


8,868 


10,805 


11,531 


10,420 


9,461 


10,591 


9,805 


9,909 


81,390 


Secondary peak height 


8,169 


10,433 


10,994 


9,938 


9,290 


10,209 


9,619 


9,682 


78,334 




Further candidate elimination stages 










< 90% of points in peaks 


8,141 


10,403 


10,968 


9,911 


9,274 


10,185 


9,605 


9,654 


78,141 


5 g' or i' points 




















Colour evolution 




















Rate of decline 




















Colour-magnitude criteria 




















Final candidates 


8,141 


10,403 


10,968 


9,938 


9,274 


10,185 


9,605 


9,654 


78,141 



light, we first define a region of the galaxy within which either the 
bulge or the disk light could be unambiguously defined. In the outer 
regions of a spiral galaxy such as M31, the visible light arises al- 
most completely from the disk. So it is possible to model the disk 
in these regions and extend the model to the inner regions of the 
galaxy. 

To greatly simplify the geometry of the galactic disk, we make 
the assumption that it is thin with an inclination of 77° (de Vau- 
couleurs 1958) and that the flux distribution is smooth across the 
disk. Another simplification we make is to essentially collapse the 
disk into a one-dimensional system. Each position within the disk 
is transformed to the semi-major axis of the ellipse that passed 
through that point. As the light from a galactic disk can often be 
modelled using a simple exponential law (Freeman 1970), the disk 
flux at any point within M3 1 can then be defined as 

Mad) = f d e- a « /a « (13) 

where fd(a,d) is the disk flux at a position within the disk with 
semi-major axis a^. As it is not possible to unambiguously separate 
disk light from bulge light, the fit is performed to the total flux data 
for ad > 40 arcmin in order to have minimal contamination from 
the bulge light. The best-fit values found for the two parameters are, 
/° = 3, 676 adu/pixel and a° d — 43.1 arcmin. 

In order to attempt to fit the bulge flux, we extend the disk 
model across the whole galaxy and then subtract the modelled disk 
light from the galactic light to leave just the bulge light and the 
residuals to the disk model. We model the bulge with elliptical 
isophotes with an axis ratio b/a — 0.6 (Ciardullo et al. 1987). We 
again transform the spatial positions of each point within the bulge 
to the semi-major axis of the ellipse on which that point lies. We 
produced the following bulge model using a standard r 1/4 law (de 
Vaucouleurs 1948, 1953): 

log[f b (a b )/fS] = -3.33[(a t /^) 1/4 - 1] (14) 



where fb{a,d) is the bulge flux at a position within the bulge with 
semi-major axis a b . A fit is performed to the bulge flux data for 
a b < 15 arcmin, so that the fit is not influenced by the disk fit 
residuals. The best-fit values are /° = 6, 914 adu/pixel and a" = 
5.1 arcmin. 



6.2 Testing the distributions 

By employing the detection probability function (Equation 12) we 
are able to test three special cases of the CN distribution in M31, 
namely that novae follow: the overall galactic light; the bulge light 
only; or the disk light only. We use the Kolmogorov-Smirnov (K-S) 
Test in all three cases to ask whether the nova distribution detected 
is consistent with being drawn from each of the three eruption dis- 
tributions. 

For the bulge-only model we use the K-S Test to determine 
whether the flux model and the CN distribution are drawn from 
the same parent population. The upper left plot in Figure 6 shows 
the cumulative distribution of bulge detection probability with in- 
creasing disk semi-major axis (ad) compared with the cumulative 
detected CN distribution with increasing disk semi-major axis. The 
K-S Test produces a probability that the two distributions are drawn 
from the same parent population of 0.43. Hence, it is clear that the 
bulge alone can give rise to the observed distribution of CNe. 

Next we test the disk-only model. The upper right plot in Fig- 
ure 6 shows the cumulative distribution of disk detection probabil- 
ity with increasing disk semi-major axis (a d ) compared with the cu- 
mulative detected CN distribution with increasing disk semi-major 
axis. The probability that these two distributions are drawn from 
the same population is 4.2 x 10~ 6 . It is therefore quite clear that 
the disk alone can not account for the observed distribution of CNe. 

Finally we test the galactic light model. The lower left plot 
in Figure 6 shows the cumulative distribution of nova detection 
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Figure 5. A plot of M3 1 flux against disk-semi-major axis distance. The 
black squares represent the flux contained within each 1 arcmin cell of the 
POINT-AGAPE data, plotted with their equivalent disk position. The grey 
points represent the M31 model flux evaluated for each POINT- AGAPE 
cell, again plotted with their disk position. The contribution to the overall 
galactic light from M32 can be seen at around a d = 70 arcmin and some 
structure - mainly from the dust lanes - can be seen within the disk across 
most of the outer galaxy. 



probability with increasing disk semi-major axis (o<j), again com- 
pared with the cumulative detected novae distribution with increas- 
ing disk position. There is a probability of 1.2 x 1CT 3 that these 
distributions are drawn from the same population. 

From the results of these tests it is clear that the bulge-only 
model does a good job of reproducing the observed CN distribu- 
tion, whilst the disk-only and galactic light models do poor jobs. 
However, given that the bulge-model over-estimates the nova dis- 
tribution near the centre of M3 1 , where the bulge dominates (see 
Figure 6 upper left plot), it is clear that a combination of both bulge 
and disk populations, with different weightings, is required to ade- 
quately model the detected CN distribution. 



6.3 The two population model 

Following the results of the testing of the three special case CN 
eruption models, it seems clear that the favoured model comprises 
a combination of both disk and bulge populations, with each pop- 
ulation having a different eruption rate per unit r' flux. To test this 
new model, we first make the assumption that the nova eruption 
probability in the disk or the bulge is proportional to the disk or 
bulge luminosity respectively: 

rb 



Pi oc a d 



ft 



(15) 



where the disk flux, / , and the bulge flux, / , are defined in Equa- 
tions 13 and 14 respectively and a d and a b are the number of CN 
eruptions per unit time per unit r' flux for the disk and bulge popu- 
lations respectively. 

In order to test this new hypothesis we define the probability 
of detecting a CN in a particular cell, given that a CN both erupts 
and is detected within one of the two POINT- AGAPE fields: 



Pi = 



(Oft + fi)- Si 

•ft 



ft 



(16) 



where 9 is the ratio of disk and bulge population eruption rates per 
unit r' flux. 

In order to constrain the favoured value of 9 we employ a max- 
imum likelihood test. The likelihood function chosen is shown in 
Equation 17 below. This function is derived from a simple Poisson 
analysis of our nova detection model, for a given underlying mean 
number of expected detections, evaluated over all possible under- 
lying means. 



r oo N 



"bin, 

n 



A"' 



(17) 



where fx is the underlying mean number of expected detections, N 
is the total number of CN detected by the POINT- AGAPE survey 
(20), rii is the number of CN detected in each data bin and Ai is the 
expected number of CN detected in each bin, given by: 



{eft + ft) ■ e* 



ft 



(18) 



In order to confine the range over which the likelihood func- 
tion is investigated we change variables from 9 to the bulge fraction 
(<E>), where <E> is defined as the fraction of the eruption probability 
within the POINT- AGAPE field due to the M31 bulge: 



4> 



2^7 = 1 



f 



(19) 



By evaluating the distribution of the likelihood function we 
derive the most likely value of $ = 0.67 and, by assuming a lin- 
ear prior in $, we evaluate confidence limits about the most likely 
value. As such, we find that the 95% confidence interval of <E> is 



0.46 sJ $ < 0.82 



(20) 



Using Equation 19, this equates to a favoured value of 6* = 0.18 
with the 95% confidence interval bounded by 9 = 0.91 and 9 = 
0.02. This result also allows us to rule out models with a d ^ er& at 
the 95% level, lending strong support to the existence of separate 
bulge and disk CN populations. 



7 M31 CN RATE 

The CN eruption probability model is used, in conjunction with the 
detection completeness data, to compute an estimate of the global 
nova rate in M31. The total number of CNe observed within the 
POINT- AGAPE fields must be proportional to the total probability 
of detecting a CN within those fields, so we have 



(21) 



where e j is the probability of detecting an erupted CN at a particular 
location, ^ is the probability of a CN erupting at that location, n 
is the number of novae detected within the POINT- AGAPE dataset 
(20), and £ is an unknown constant relating the detection probabil- 
ity to the recovered number of novae. The definition of the eruption 
probability given in Equation 15 is used to define 



eft + ft 



fb 
■I 3 



(22) 
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Figure 6. The major axis distribution of the 20 detected POINT-AGAPE novae (grey histograms) along with theoretical predictions of four models (black 
lines). The upper left panel shows the bulge-only model, whilst the upper right panel shows the disk-only model. The lower left panel shows the galactic light 
model, whereas the lower right panel represents the most probable distribution (see Section 6.3). 



The value of the multiplier £ can be computed for a number of 
different values of 6, thus producing a range of M31 nova rates. 
However, we will restrict the values of 9 examined to those that re- 
late to specific physical situations: 9 — 0, the bulge-only system; 
6 = 1, the galactic light scenario; 9 — ► oo, the disk-only system 
and 9 = 0.18, the favoured value produced by the maximum like- 
lihood analysis of the two population model. The computed values 
of £ for these models are given in Table 3. 

The global M3 1 CN number can now be computed from 

N = - (23) 

<P 

where N is the global nova number and tp is a constant multiplier 
that accounts for the proportion of the total galactic eruption prob- 



ability that has been sampled by the POINT-AGAPE survey, tp is 
defined by 

In order to evaluate the sum over the entire galaxy (the denominator 
of Equation 24) we extend the 1 arcmin grid, initially used to model 
the completeness, over all space. As galactic disks are known to 
be truncated radially at a distance of 3^1 scale lengths (van der 
Kruit & Searle 1981; Pohlen et al. 2000), we evaluated the sum 
out to a distance of 20 kpc (Ibata et al. 2005) (equivalent to a disk 
semi-major axis distance of ~ 90 arcmin) from the centre of M31. 
Table 3 shows the computed values of tp for the models tested. We 
are now able to compute the global number of CNe that erupted in 
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Table 3. The computed values of the nova rate normalisation (£), the approximate bulge-to-disk probability ratio, the fraction of novae erupting within the 
POINT-AGAPE fields (tp), the underlying number of nova eruption (AT) during the survey lifetime, the M31 bulge nova rate (Nbulge), the disk nova rate 
(Ndisk) an d the global nova rate (N) for a range of different CN eruption probability models. 
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0.31 


233 




82 ± 18 


82 ± 18 



M3 1 during the POINT- AGAPE observing baseline. These values 
are also shown in Table 3 for the four model examples. 

However, to calculate the global M3 1 CN rate, we first need 
to take account of the finite observable lifespan (as defined solely 
by our observations) of each of the detected novae. As we did not 
require that a nova's light-curve should be completely contained 
within our data, our effective baseline for each CN is extended by 
the lifespan of that particular CN. As the novae have been seeded 
uniformly over the POINT-AGAPE fields and the seeded nova is 
selected randomly, the baseline of observations has been extended, 
on average, by the mean lifetime of all 20 novae: 

T — Tfjaseline tnova (25) 

where, TbasMne is the time between the first and last POINT- 
AGAPE observation (2.472 years), and i„ ova is the mean lifetime 
of the 20 POINT-AGAPE novae (0.359 years). The effective base- 
line for CNe of the POINT-AGAPE survey is therefore 2.830 years. 
The nova rate for M3 1 is then 

N=4f- (26) 

The computed M3 1 global CN rates for our four model scenarios 
are given in Table 3. This illustrates that the predicted overall nova 
rate is only modestly dependent upon the eruption model chosen. 
However, the separate bulge and disk novae rates show a strong 
dependence upon 8 as expected. 

The errors shown for the nova rates in Table 3 are generated 
solely from the Poisson errors related to the size of our nova cat- 
alogue. Given the small size of this catalogue, this source of error 
(~ 22%) is expected to dominate over all others. The other main 
sources of error arise from the completeness calculations, the lack 
of fast novae in the catalogue, the possible misidentification of no- 
vae and the modelling of the surface brightness. In addition, our 
limited knowledge of the internal extinction of M31 makes it diffi- 
cult to estimate the errors introduced into the completeness by its 
exclusion. However, we expect these errors to be small. The max- 
imum r' extinction expected in the disk is ~ 0.7 magnitudes; as 
all the POINT-AGAPE novae were followed through at least one 
magnitude, it is expected that few, if any, novae were missed due 
to extinction problems. Should a non-CN have been wrongly in- 
cluded within the nova catalogue, this would directly result in a 5% 
reduction of the nova rate (see Equation 21), coupled, indirectly, 
with a further maximum decrease of 5% from the completeness 
model. Hence, the misidentification of a single nova generates a 
maximum absolute error of ~ 7%. However, the errors introduced 
into the completeness through misidentification are highly likely to 
be much less than 5%. Likewise, a single CN which is "missed" 
due to extinction would induce a maximum increase of ~ 7% of 
the global nova rate. Although there is some error in our modelling 



of the M31 surface brightness within the POINT- AGAPE fields, the 
good fit of the model shown in Figure 5 suggests that this is modest. 

Taking the above error discussion into account, we can use the 
value of 9 = 0.181^10 deduced from the likelihood analysis of the 
two population model to produce a model constrained estimate of 
the true nova rate of M3 1 . By evaluating the error on the bulge frac- 
tion determination we are able to show that the Poisson errors are 
still the dominating error source although there is also a significant 
contribution from the uncertainty in the model. Hence we arrive at 
the following estimate of the true observable nova rate of M3 1 : 

A> = 65±if year" 1 (27) 



8 DISCUSSION AND CONCLUSIONS 
8.1 Novae as distance indicators 

Sections 3 and 4 report the calibration of both the MMRD and tis 
relationships using the POINT-AGAPE CN catalogue, which in- 
cludes a recalibration of each nova's decay rate and speed class and 
assessments of the uncertainties in the maximum light and line of 
sight extinction to each nova. 

The two POINT-AGAPE MMRD relationships (see Equa- 
tions 5 and 6) are consistent with the existence of an MMRD re- 
lationship for the r' and i' filters. In fact, the observed scatter in 
both relationships, although higher than that of recent Galactic cal- 
ibrations, can be accounted for solely by extinction and maximum 
light uncertainties. We are able to show that, for the speed classes 
used for the MMRD calibration, missing the maximum light of a 
nova by up to a week is not the dominating factor in the MMRD 
scatter. We also find that for the bright novae for which a linear 
MMRD is expected, extinction corrections make little difference 
to the slope determination. However, it is clear that a better under- 
standing of the extinction affecting the POINT-AGAPE novae is 
required in order to make a more precise calibration of the MMRD 
within M3 1 . Little more can be said about the comparison between 
previous MMRD relationship calibrations and the POINT- AGAPE 
calibrations, as the Galactic MMRD (and previous M31 relations) 
are calibrated using bluer filter bands than the POINT- AGAPE fil- 
ters. In fact our calibrations constitute the first attempt to do so 
using Sloan filters. Given that CNe become bluer as they decline, 
we would naively expect the POINT-AGAPE r' and i' slopes to be 
steeper than the galactic V-band slope, whereas we find that the r' 
slope is remarkably similar, with the i' slope being much steeper, as 
expected. However, it should also be noted that the r' filter contains 
the Ha emission line. As CNe are known to remain bright in Ha 
long after the visible light-curve has diminished, this may be ad- 
versely increasing our measured r' decline times. It is also known 
that the decline of the Ha emission of a CN is not well correlated 
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with its Ha luminosity at maximum (Shafter 2005). As such, this 
could potentially detract from the usefulness of any r' MMRD re- 
lationship (Shafter, private communication). 

The MMRD relationship may be be used as a tool to measure 
the relative distance between two populations of novae. However, 
given that our calibrations are the first to be carried out for the Sloan 
r' and i' bands, it would be inappropriate to attempt to estimate the 
M3 1 distance by comparison with Galactic V and B band relation- 
ships. 

The analysis of the r' and i' tis relationships within the 
POINT-AGAPE catalogue is, like the MMRD relationship, dom- 
inated by the extinction uncertainties within the data. However, the 
extent of the scatter observed in both the r' and i! data can not be 
accounted for by the extinction and maximum light uncertainties 
alone. A comparison of our tis values with those for bluer bands 
are, however, consistent with a CN becoming bluer following max- 
imum light. 

Figure 3 shows the distribution of t n scatter, n days after max- 
imum light. This plot clearly indicates that the scatter between the 
light-curves is large over the entire period sampled. Also it is clear 
that there is no evidence of a minimum in the scatter for times 
around 15 days. The small minima at ~ 30 days in the r' data 
and ~ 35 days for the i! are related to the sampling of the surveys. 

To some extent, the tis analysis is limited by the temporal 
sampling of the POINT-AGAPE survey. Unlike the MMRD rela- 
tionship that requires good sampling around the peak of the light- 
curve and some good sampling of the subsequent decline, to test the 
tis relationship one also requires good sampling of the light-curve 
specifically at ~ 15 days after peak. Due to the make up of the sur- 
vey, the light-curves are generally constructed from short periods 
of good sampling, followed by regions with no data (see Figure 2 
of Paper 1). As such, a relatively large amount of extrapolation is 
required to estimate each nova's flux between observations. Given 
the rather erratic behaviour of a CN's light-curve, the estimation of 
the errors induced by linearly interpolating over large periods with 
no data is a far from trivial task. The r' data are again likely to be 
adversely affected by the Ha emission. We can conclude that the 
POINT- AGAPE CN catalogue shows no evidence of a tis relation- 
ship, nor strong evidence of convergence at another timescale. We 
would require nova light-curves with much more uniform sampling 
than the POINT- AGAPE novae to be able to make a more definite 
statement regarding the t n relationship's overall validity in the r' 
and i' filters and its potential usefulness. 

8.2 Completeness 

Overall, the method employed to evaluate the completeness of the 
POINT-AGAPE CN catalogue generated by the nova pipeline al- 
lowed us to obtain a very good understanding of the CN detection 
efficiency of both the survey and the pipeline. The completeness 
analysis took into account a variety of possible selection effects 
which prevent us from detecting novae. These selection effects in- 
cluded the strongly varying surface brightness of M31, the range of 
morphologies exhibited by CN light-curves and the temporal sam- 
pling of the POINT- AGAPE survey. 

Until very recently the detection of novae relied solely upon 
visual detection, often by the "blinking" of images. Even the most 
recent surveys (Shafter & Irby 2001, for example) have relied on 
some visual inspection, particularly to aid in the detection of the 
faintest novae. The majority of past nova surveys have also re- 
lied upon visual inspection of light-curves to determine the like- 
lihood that an object was a CN (e.g Ferrarese et al. 2003). As our 



pipeline uses much more robust methods and objective selection 
criteria to both detect and classify potential CNe, we are confident 
that the completeness of the catalogue is well understood. Whilst 
the POINT- AGAPE CN catalogue may not be complete 4 , we are 
nonetheless able to quantify our completeness. 

There were, however, a number of factors that have not been 
taken into account by the completeness analysis. As was discussed 
in Section 2, our knowledge of the internal extinction of M3 1 is lim- 
ited and these extinction uncertainties have not been built into the 
completeness computations. Whilst the extinction may be dimin- 
ishing our ability to detect novae, especially fainter CNe, its rela- 
tively small magnitude should not be too troublesome. There are 
no very fast novae (t2 ^ 10 days) within the POINT- AG APE cat- 
alogue . Whilst this may simply be indicative of our small sample 
size, there may also be additional selection effects - due to the tem- 
poral sampling of the POINT- AGAPE survey - that are preventing 
us from detecting novae of this class. Both of these effects poten- 
tially prevent us from observing novae erupting during the survey. 
Consequently, the computed completeness is likely to be an over- 
estimate. Given the form of the extinction within M31, it is also 
likely that the completeness has been over estimated to a greater 
extent within the disk (due to its generally greater extinction) than 
within the bulge. However, Shafter & Irby (2001) used observa- 
tions of the planetary nebula population of M3 1 to conclude that 
(in Ha) the observed CN population is not significantly affected by 
extinction. As it is expected that fast novae are more likely to be 
observed within the disk of a galaxy (see Section 1.4), the probable 
exclusion of very fast novae from the catalogue again leads us to 
conclude that, if anything, we have over-estimated the complete- 
ness and so under-estimated the nova rate. 



8.3 M31's CNe population 

The analysis of the observed CN distribution within M3 1 allows us 
to develop a basic model of the underlying CN distribution. Fol- 
lowing the separation of the disk and bulge through modelling of 
the surface brightness, we are able to show that, within the two 
POINT-AGAPE fields, the observed CN distribution does not fol- 
low that of the galactic light. Nor is the observed distribution likely 
to arise solely from the disk. Although the bulge alone can support 
the observed distribution, a combination of a bulge and disk popu- 
lation seems required to fully reproduce the observed distribution. 

A maximum likelihood analysis of the two population model 
(see Section 6.3) indicates that the ratio of the disk and bulge pop- 
ulation eruption rates per unit r' flux is 0.18. This result is consis- 
tent with previous findings (Ciardullo et al. 1987; Capaccioli et al. 
1989; Shafter & Irby 2001) which reported that the M31 novae are 
primarily associated with the bulge. Shafter & Irby also reported an 
eruption rate per unit B flux within the bulge of up to an order of 
magnitude greater than that of the disk. Figure 7 shows a schematic 
plot of our "best" CN eruption model. 

The range of bulge-to-disk eruption rate ratios that are consis- 
tent with the range of observed CN distributions (0 = 0.18tg;?o) 
leads to a range of expected bulge-to-disk nova eruption rate ratios. 

4 A number of CN candidates contained within the POINT- AGAPE dataset 
are known not to be contained within our catalogue (An et al. 2004; Feeney 
et al. 2005) for reasons we understand (Feeney et al. 2005). 

5 PACN-99-07 has a great uncertainty in its speed class assignment. Al- 
though initially classified as a very fast nova, it is thought more likely to be 
a moderately fast or slow nova. 
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Figure 7. The favoured M31 eruption probability model (9 = 0.18), over 
plotted with the position of the POINT-AGAPE fields and the 20 detected 
CNe. 



The global bulge-to-disk CN ratios range from 5 : 1 for a bulge 
dominant population (8 — 0.1) through to 1 : 4 for the distribution 
following the surface brightness (6 = 1). The expected nova ra- 
tios within the POINT- AGAPE fields themselves range from 20 : 1 
(0 = 0.1) to 2 : 1 (0 = 1). However, the POINT- AGAPE survey 
of M3 1 covers a much greater surface area of M3 1 than all previ- 
ous nova surveys (see Shafter 2005, for a summary), which have 
concentrated mainly on the bulge. The POINT-AGAPE survey has 
given us much better coverage of the M31 disk and its CN pop- 
ulation. As a result, these previous surveys will have all observed 
a distribution that appears to be much more bulge dominated than 
that of the POINT- AGAPE catalogue. 

The analysis of the M31 CN distribution is, however, limited 
by a number of considerations. These are mainly the small size 
of the POINT-AGAPE CN catalogue, the simplistic nature of the 
M31 surface brightness models and the uncertainties arising from 
the completeness modelling. Further, the CN distribution analysis 
relies upon the modelling of the M31 surface brightness. However, 
this modelling only includes the "normal" disk component and the 
bulge component of the surface brightness. Other components, such 
as the spiral arm structure, the dust lanes (and extinction within 
M31 in general) and M32, are not taken into account. As with the 
completeness analysis, these effects are expected to have greater 
affect within the disk than the bulge. Thus, the inclusion of the 
dust lanes and spiral structure within the models could lead to an 
increase in the expected number of disk novae. 



8.4 The M31 and Milky Way nova rates 

By extending the M3 1 CN eruption model over the whole galaxy 
we are able to produce an estimate of the global nova rate. The com- 
puted global observable nova rate of M31 is 65+ ^ year -1 , with a 
bulge rate of 38±\l year -1 and a disk rate of 27±{ 9 5 year -1 . This 
result is at the limit of being consistent with that of the most robust 



previous calibration, which found a global rate of 37t 8 year -1 
(Shafter & Irby 2001). However, our results are much higher than 
all previous results, including the Shafter & Irby determination and 
those of Hubble (1929) (~ 30 year -1 ), Arp (1956) (24±4 year -1 ) 
and Capaccioli et al. (1989) (29 ± 4 year - ). The ratio between the 
bulge and disk nova rate is also markedly lower than that computed 
by Shafter & Irby's (from their maximum likelihood analysis of the 
M31 CN distribution). 

Despite its apparent high value, we are confident that our com- 
puted nova rate is a true evaluation of the nova production rate of 
M31. Our robust completeness analysis and objective selection cri- 
teria lead us to believe that the completeness of previous surveys 
may have been over estimated. Also, given their bulge-centric na- 
ture, many previous surveys are likely to have under estimated the 
contribution from disk novae. Sources of concern in our estimated 
rate arise again from the extinction uncertainties and the lack of 
very fast novae in the catalogue. However, both these factors po- 
tentially lead to a further increase in the predicted rate, so we are 
led to conclude that the true global CNe rate of M3 1 is higher than 
was previously thought by around 50%. 

We can use our estimated global M3 1 nova rate, along with our 
computed eruption rates for the bulge and the disk, to produce an 
estimate of the global Galactic nova rate. Using a similar method to 
that outlined in Shafter (2002), we assume a Galactic disk-to-bulge 
luminosity ratio of ~ 8 and a Milky Way to M3 1 luminosity ratio 
of 2/3. Hence, we compute a global Galactic nova rate of 34^{l 
year -1 , with a disk rate of 201" year -1 and a bulge rate of 14±g 
year -1 . These rates are broadly consistent with other Galactic es- 
timates based on M3 1 data and broadly consistent with estimates 
based upon direct observations of Galactic novae (Shafter 2002). 
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